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We consider overdamped diffusion processes driven out of tfiermal equilibrium and we analyze 
tiieir dynamical steady fluctuations. We discuss the thermodynamic interpretation of the joint 
fluctuations of occupation times and currents; they incorporate respectively the time-symmetric 
and the time-antisymmetric sector of the fluctuations. We highlight the canonical structure of the 
joint fluctuations. The novel concept of traffic complements the entropy production for the study 
of the occupation statistics. We explain how the occupation and current fluctuations get mutually 

' coupled out of equilibrium. Their decoupling close-to-equilibrium explains the validity of entropy 

' production principles. 

o ■ 

(N 
> 

^ . I. INTRODUCTION 

0^ ■ Thermodynamics already ceases to be exact under the microscope. That line of thought appeared a century ago 
explicitly so in the introduction to Einstein's paper on the atomic hypothesis, The ensuing theory of Brownian 
' ' . motion started a dynamical fluctuation theory, and essentially the same mathematical model was further used by 
O ' Onsager and Machlup, to describe small macroscopic fluctuations around the relaxation to equilibrium. It has 
^ , developed into the (sometimes called Lagrangian) approach to noncquilibrium statistical mechanics which is based on 
giving a weight to available trajectories in terms of an action that integrates a space-time local Lagrangian. Since then 
various methods have been developed to extract from that action functional physically relevant information about the 
d • steady state statistics. 

^ [ In recent years there has been a revival of noncquilibrium fluctuation theory, with many contributions on various 
. aspects of the problem. One class are static fluctuations referring to the statistics of macroscopic fluctuations upon 
' a nonequilibrium stationary distribution; they have been mostly studied for driven lattice gases in the hydrodynamic 
^ , regime, 0] , by generalizing the Hamilton- Jacobi method originally proposed in [J . Related yet fundamentally different 
are the dynamical fluctuations concerning time-integrated observables. That is what this paper is about. Wc will 
introduce and also compare our present work with other recent approaches in Section lll CI There has been a particular 
interest in the (integrated) current and entropy production fluctuations, also in the context of the celebrated fluctuation 
^ [ symmetry. It is nearly impossible to mention all the relevant references; as an example we include [ Si, A Ji, .8., j),, 10, Ijj. 

An extension of the dynamical approach to the time-symmetric domain was initiated in [13, particular, 
' it has been shown how the nonequilibrium fluctuations of the steady state occupation times are related to the mean 
^ , entropy production, [l^ . That relation proves useful for understanding the status and the limitations of entropy 
' production principles in characterizing the steady state, [l5l|. As further unfolded in [l6l | within the framework 
OO . of Markov jump processes, a particularly simple and generic structure of dynamical fluctuations emerges when both 
' the time-symmetric (occupation times) and the time-antisymmetric (current) fluctuations are observed jointly. Their 
I correlation can be neglected in a close to equilibrium regime, where the familiar linear irreversible scheme based on the 
. entropy production alone exists. However, the coupling between the time-symmetric and time-antisymmetric sectors 
' becomes essential far from equilibrium. That can be seen as a fundamental reason for highly complex patterns of 
. e.g. current fluctuations which arise when observing their marginal distribution alone. The main question remains 
' to see some systematics in the fluctuation formula; and to identify terms in the rate functions that have a general 
thermodynamic meaning. In order to complete that program the following examples will help as orientation. 

The plan of the paper is as follows. In Section[lllwc introduce a class of driven diffusion models on which we formulate 
our main questions and answers. In Section IlIII wc define the dynamical entropy, composed of a time-antisymmetric 
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■ part (entropy flux) and a time-symmetric part (dynamical activity or traffic). Our main results come in Section HVl 
where we show how entropy production functionals govern the dynamical fluctuations of diffusion processes. We see 
then how the time-symmetric and the time-antisymmetric fluctuations get coupled, and we discuss the validity of 
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some entropy production principles. The mathematics involved is the standard Ito/Stratonovich stochastic calculus 
complemented with the path-integral formalism; for convenience we give a brief review in the appendix together with 
some additional references. 



II. MODEL AND MAIN RESULTS 

We start with an example of a one-dimensional overdamped Langevin dynamics driven from equilibrium by a 
nongradient force. This is further generalized to a larger class of diffusions, and in arbitrary dimension. 

A. Basic example: diffusion on a circle 

We consider a particle undergoing an overdamped motion on the circle with unit length. This means that we ignore 
inertial effects, so that forces will be proportional to velocities rather than to acceleration. The particle moves under 
the influence of a stochastic force (noise) , because we imagine the system to be connected to an environment at inverse 
temperature p. Added to that, there are also deterministic forces: a periodic potential landscape U{x), and a periodic 
(nongradient) force F{x), which will drive the system out of equilibrium. Note that on the circle every function F 
is periodic, but is not in general a derivative of a (periodic) function (for example if > everywhere). With some 
(possibly inhomogeneous) mobility xi^) > this is modeled by the Ito-stochastic dynamics (see Appendix A): 



dxt - xixt)[F{xt) - U'ixt)]dt + D'{xt)dt + ^2D{xt)dBt (1) 

where dBt is standard Gaussian white noise, so that the process is Markov. The prime as superscript is a shorthand 
for the spatial derivative. The diffusion coefficient is D = xIP: in agreement with the condition of local detailed 
balance. The drift counter-term proportional to D' then ensures that the case = is an equilibrium dynamics, see 
below. For the following we assume that this stochastic process ([T]) always relaxes to a unique stationary distribution 
which mathematically amounts here to smoothness conditions on _F, U and x, see e.g. in |17| . 
The corresponding Fokker-Planck equation for the time-dependent probability density /ij is 

^ + j;. (^) = 0, J, = x^^{F ~U')- (2) 

where is the probability current as the sum of a drift and a diffusion component. 

The stationary density, called p in the sequel, solves the stationary equation j'^ = 0, i.e., 

XPiF - U') - Dp' = jp (3) 
is a constant. For F = 0, equation Q has the solution 

p(.T) = ie-^^(-), Z= l\-P^dx (4) 
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and the corresponding stationary current is jp = 0; this is a detailed balanced dynamics with p the equilibrium density. 
When adding a nongradient driving force, Fdx =^ 0, the stationary density obtains the form 

where 

F dz for y < X 

jy Fdz + Jq Fdz for y>x 

is the work performed by the applied forces along the positively oriented path y x. In this model the stationary 
current can be computed by dividing the stationary equation ^ by px and by integration over the circle: 

'''' /o(PX)"^da; 

where W = Fdx is the work carried over a completed cycle. The non-zero value of this stationary current indicates 
that time-reversibility is broken. In the simplest nonequilibrium setting when U = Q and F,x > are some constants, 
the steady state has the uniform density p{x) ~ 1 and the current is jp = xF. 



B. General model 



As a generalization of the driven diffusion on a circle, we consider a class of d— dimensional inhomogeneous diffusions 
introduced by the following equation, which has to be interpreted in the Ito way (see Appendix A): 

dxt = {xixt) [Fixt) - VU{xt)] + V ■ Dixt)}dt + y/2D{xt)dBt (8) 

The mobility x(^) is now a strictly positive (symmetric) d x d— matrix which is related to the diffusion matrix D{x) 
as xi^) = (3D{x) with some fixed homogeneous inverse temperature /3 > 0, and the d— dimensional vector dBt has 
independent standard Gaussian white noise components. The notation V • D stands for the vector with components 
dryDya, i-c, with the derivative acting on the left indices of a matrix. From now on we also use the dot to denote 
scalar product, e.g., F-G = F-^G-^, whereas no dot is used for a matrix acting on a vector, e.g., (x F)a = Xa-yF^. 
We restrict ourselves to two types of boundary conditions: 

(1) periodic — the particle moves on the unit torus [0, 1)'' and the fields U , F, and x ^re smooth functions on the torus; 

(2) decay at infinity — the potential U grows fast enough at infinity so that the particle is essentially confined to a 
bounded region, i.e., the density and its derivative vanish at infinity. 

Under cither of the above boundary conditions wc can simply ignore boundary terms when performing integrations 
by parts. The particles are essentially confined in their configuration space. 

The probability density evolves according to the Fokker-Planck equation associated with ([5]) : 

^ + v-jV,, -0, J^ = x^^iF-vu)~DV^l (9) 

where again x = PF>. The stationary condition reads V ■ jp = 0; in contrast to the previous one-dimensional example 
the stationary current in general becomes inhomogeneous here. Equation Q can be interpreted as giving the evolution 
of the density profile /it of a macroscopic amount of independent particles each moving according to ([S]): at times 
t > and for all observables /, 



f{xt)) = / f{x)iit{x)dx (10) 
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where is the average over the path-space distribution started at hq. 

Moreover, the then also gives the expected profile of the 'real' particle current at given density fi: again under the 
diffusion process started from density /zq, 

f{xt)odxt)^^^ dt j f{x)j^^{x)dx (11) 

which involves the Stratonovich-stochastic integral; see (|A8p in Appendix [X] for its derivation. Taking formally 
/( - ) = 5( ■ — x), we indeed recognize jit and j^^ as the time-dependent (or transient) local density, respectively the 
local current for that density. 



C. Questions and first answers 

Since the Langcvin dynamics ([1]) or ([8]) is assumed to have a unique stationary distribution, the stationary density 
p can be measured as the most typical time-average of the occupation over a large time interval. Any different 
measurement outcome, specified by some density /i, is well possible but it is a rare event with an increasingly small 
probability to be observed as the duration T f +oo. The large deviation theory, [3, [l^, tells us that the generic 
asymptotic law for such fluctuations is an exponential decay, PxifJ^) ~ g-'^^M with some rate function /(/i). Here 
comes a first question: 

QA: How can we find the rate function I{p.) and what is its thermodynamic meaning? 

It turns out that I{p.) can be expressed in terms of the traffic measuring the dynamical activity in the system, and 
being introduced as a time-symmetric counterpart to the entropy production. We also give a relation between the 
traffic and the entropy production which appears to be general for diffusion processes. 

A very similar question can also be asked about the current statistics, with rate function but the study of /(/i) 
is quite new. 



QB: What arc the fluctuations of the time- averaged current around its most probable value jp, and how are they 
correlated with the fluctuations of the occupation times? In particular, what arc the joint fluctuations for both 
time-symmetric and time-antisymmetric observables? 

To answer these questions, we give an explicit form of the dynamical fluctuation functional such that 

-F'(Mjj) ~ e"^'^^'''-'-' is the joint probability to see both /i as the statistics of occupations times and j as the 
time-averaged current. The functional I{p,,j) has a general canonical structure and it is given in terms of the entropy 
production and the traffic. Starting from this basic functional other fluctuation formulas can be obtained by the 
contraction principle familiar from large deviation theory. For example, one recovers I{fi) by minimizing I{p,j) over 
all admissible currents j. 

The above questions and the methods that we take below are not entirely original. They have appeared in the 
mathematical literature in a systematic way since the theory of large deviations was introduced in the framework 
of Markov diffusion processes, see [l^, [l^, The relevance to physics and to statistical mechanics in particular 
is obvious, but the thermodynamic interpretation of the resulting dynamical fluctuation functionals has not been 
systematically investigated. A first study can be found in Our paper will add to that, starting from the next 
section. There have of course been many other studies of dynamical fluctuation theory in the literature. We mention 
in particular the works of Derrida and Bodineau, see d, Q and of Bertini et al, e.g. in 1,0). As we are dealing here 
with diffusion processes, our approach is especially similar to what one is doing for a macrostatistical theory, where 
the hydrodynamic fluctuations can be viewed as solution of some infinite-dimensional diffusion process. Possible 
differences with the existing work are first of all that the problems related to the diffusion-approximation or to a 
hydrodynamic rescaling do not enter in our work. We just start from a finite-dimensional diffusion process as such, in 
the same way as we could start from a Markov jump process, boundary or bulk driven, and without extra rescaling. 
In other words, we prefer to split the problem of hydrodynamical scaling with possible diffusion approximation from 
the problem of studying the dynamical fiuctuations. Secondly, we are concerned here directly with the stationary 
fluctuations, without passing via the transient regime. Thirdly, as soon will become clear, we like to emphasize a 
(thermodynamic) canonical structure in the joint fluctuations of density and current. In order to achieve that, it is 
natural to introduce the novel concept of traffic (as in the next section) whose thermodynamic interpretation remains 
to be fully elucidated. Our three main results are in Sections IIVBI IIVDI and IV Al 

III. THERMODYNAMIC CONSIDERATIONS 

In this section we introduce the two main players in the presented fluctuation theory: the entropy production and 
the traffic. While the former notion is rather standard, the latter appears to be new. 

A. Entropy production 

The entropy production is the change of entropy in the world between two moments in a process. The very notion 
of entropy is however controversial when trying to extend it to noncquilibrium situations. The problem is mainly that 
we do not have a good understanding of the relevant quantities or of the macroscopic variables that characterize the 
noncquilibrium state and its evolution. Nevertheless we can try to work with the notion of entropy production as it 
has come to us from considerations of heterogeneous or of local equilibrium. There are two basic contributions to 
that entropy production. 

The entropy of a closed macroscopic system is defined via counting microstates, which at the same time gives both 
the static fluctuation theory with the entropy as a fluctuation functional and a second law inequality, |2l| . For a 
mesoscopic system this cannot be applied as such. However, one can think of an ensemble of A'^ — > oo independent 
copies. A density fj, then becomes a macro-observable telling us the relative occupation of the space, and the associated 
counting entropy equals the relative entropy (with respect to the flat distribution): 

s(/x) = — y ^j.{x) log fj,{x) dx (12) 

By construction it is the (static) fluctuation functional in the probability law for observing the empirical density fi 
when sampling the particles from the flat distribution. For F = \7U = (i.e. for a thermodynamically closed system) 
the entropy rate satisfies ds/dt > along the solution of the Fokkcr-Planck equation. 

An open system dissipates heat that results in an extra entropy production in the environment. The rate Q(/i) of 
mean heat dissipation comes from the work performed by the nongradient force and from the change of the energy of 



the system: 



= j{F-VU)- dx + Jv- (Uj^) dx 



(13) 



where we have used the Fokker-Planck equation ([9]) . Under either of the two boundary conditions considered in this 
paper (i.e. periodic or decaying at infinity, see Section fll Bp . the last integral equals zero. That argument will often 
come back in what follows. 

The environment is a heat reservoir at inverse temperature (3 that remains itself at equilibrium during the whole 
process (the weak coupling assumption). Hence, the mean entropy flux equals (3Q{fi) and the total entropy production 
rate reads 

a(M)=/3Q(/i) + ^(/.) 

= JiPF- PVU -V\og^i) ■ dx (14) 

jf, ■ {fioy^jf, dx 



where we have used again the Fokker-Planck equation @ and the boundary conditions. In the language of irreversible 
thermodynamics, (3F — \7\og{fie^^) is a generalized thermodynamic force conjugated to the current j^. Clearly, 
'^{p) ^ in agreement with the second law, proving the thermodynamic consistency of our diffusion model. 

B. Dynamical entropy and traffic 

The diffusion as defined by the stochastic equation ([H]) is a Markov process. Its randomness can be characterized 
by the dynamical entropy which is a dynamical variant of the static entropy p2p obtained by replacing the density 
/i with the path-space distribution of the process. Such a construction is often useful in the theory of dynamical and 
stochastic systems; here we review how it is linked to the thermodynamic entropy production defined in the previous 
section and what other information it provides. 

We start by writing the density of the path-distribution of the process with respect to a suitable reference 
process P", for which we take the one corresponding to F = VU ~ (see Appendix B): 

dP^„(a;) =/.to(a;o)e-^"''^^^"''''*"^«^^"^""'^"'MP°(w) (15) 

where w = {xt)J^Q is a trajectory and /zq is an initial density. The weight in the exponential is split in a time-symmetric 
part 

C+{x) = ^[(F - VC/) • D{F - VC/)] + ■ [D{F - VC/)] (16) 

and a time-antisymmetric part 

C~{x,dx) = -^{F -VU)odx (17) 

The dynamical entropy is the relative entropy of the path-distribution P^ with respect to the reference P*^, over 
some fixed time interval [0,T]: 



^M = (log^) =-^(mo)-(/ C+{xt)dt) -(/ C-{xtAxt)) (18) 



The time-antisymmetric contribution to the dynamical entropy equals, applying the identity (jlip . 



^ C-{xt,dxt))^^ = ^ dtJ{F-VU)j^,dx^!^J^ Q{^it)dt 



(19) 



Hence, it is given in terms of the entropy flux. This is an instance of a fairly general observation that lies behind a 

famous fluctuation symmetry of the entropy production, [2^ . [23l . [23 | . 

The time-symmetric contribution is computed analogously via (jlOp to obtain 

( / C+{xt)dt) = [ T{pit)dt (20) 
where we have introduced a new quantity T{fj,), called traffic^ which equals 

"^(m) = y / " ■ ■^^^ " + f / ' '^^^ " ^^^^ ^^^^ 

By construction, the traffic is that part of the dynamical entropy which originates from the time-symmetric fluctua- 
tions. Finally, we arrive at the next two equivalent expressions for the dynamical entropy 1)18^ : 
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where we have inserted (HH). 



I?(Mo) = -s(Aio) - / T{fit)dt+^ [ Q{fit)dt 
Jo ^ Jo 

= -^[s(a'o) + s(mt) + 2 / T{nt)dt- f a{fit)dt] 



(22) 



C. Traffic versus entropy production 



A remarkable and important feature of diffusion systems, not valid in general beyond the diffusion approximation, 
p^ . is that the traffic (PT]) and the entropy production (fH|) are not independent of each other. Indeed, one checks 
the relation 

Ti^i) = ^ + J (23) 

in which the last term only depends on the distribution fi and not on the imposed potential U nor on the driving 
F. As we will see in the next section, differences in the traffic and in the entropy production when varying U and 
F determine the asymptotics of dynamical fluctuations. Hence, relation ()23p brings about a simplification in the 
structure of fluctuations that is characteristic and restricted to diffusions. We believe that this also indicates that a 
more general nonequilibrium theory should also reach beyond the Langevin or diffusion approximation. 



IV. JOINT OCCUPATION-CURRENT STATISTICS 



According to Einstein's fluctuation theory, the entropy and derived quantities govern the structure of static fluc- 
tuations. In this section we explain how the time-symmetric and time-antisymmetric components of the dynamical 
entropy play an essential role in dynamical fluctuation theory. 



A. Definitions 



A basic and time-symmetric dynamical observable is the empirical density of the occupation times defined as the 
fraction of time spent in the neighborhood of every point x, over a fixed time interval T: 

Mx) = ^ / 5{xt-x)dt (24) 
^ Jo 

This is a path-dependent observable as it varies over the paths lu = (xt)^Q. For the large time asymptotics T — > oo, 
we have ^lt —> p almost surely and independently of the initial condition. We will be concerned with physically 
determining the fluctuation functional /(/i) that enters the large deviation law 



P(/iT = m) - e-"('^) 



(25) 



That has to be understood as an asymptotie formula that becomes an equahty after taking the logarithm and 
dividing by T ^ cxd in both left- and right-hand sides; see p^ . [25l | for a more precise mathematical formulation. 

The time-antisymmetric observable of special relevance is the empirical current, formally defined as 

1 

3t{x) = 7^ I 5{xt~x)o dxt (26) 
^ Jo 

It depends again on the (random) path oj = {xt)f^Q and it measures the time-averaged current while in x, as in 
Its steady average equals jp] moreover, jt — > jp for T oo almost surely. The rate function I{j) is defined 
analogously to (PS)) . but for the probability P{jT = j)- 

Both /(/-i) and /(j) can be obtained in principle from the joint fluctuations. The large time asymptotics of the joint 
fluctuations of ftx and jx is described by the fluctuation functional I{fi,j) such that 

P(/2t = m; Jt = j) e-"('^^J) (27) 

always logarithmically when T — > oo. A first observation is that = oo = I{j) whenever j is not stationary, i.e. 

for V ■ j 7^ 0. Indeed, for any smooth bounded function Y one has 

jYV-jrdx^-^J^ VYixt)odxt^-^[Y{xT)-Yixo)]^0 (28) 

and hence, in a distributional sense, V • jt — > for T — > oo along any particle trajectory, which proves the above 
statement. That is why from now on we always assume that V • j = 0, unless otherwise specified. 



B. Result 1: traffic and entropy production determining the joint fluctuations 

To compute we use a standard large deviation technique, sometimes referred to as Cramer tilting, see e.g. 

[igl . [20t . We modify the driving of the original dynamics changing F into some G and we take care that this new and 
modified dynamics is chosen so that fj. and j become both stationary (and hence are typical as large time averages). 
The potential U remains untouched. The more formal argument goes as follows. 

Let Pp and Pq be the stationary path-space distributions for driving F respectively G, both processes started at 
their respective stationary distribution. We refer to Appendix |B] for that notion of path-space distribution. We must 
compute the probability -P_f[A] of the event A containing all the paths ui over time T such that iIt = and jt = j, 
see (HZl): 

dPaii^) (^) ^['^^ = JT = j] (29) 

The density of path-space distribution Pp with respect to path-space distribution Pq can be written out explicitly 
starting from (|B3p : we find 



where we have made use of definitions and pT|) . The point is now that this can be fully expressed as a function 
of the path-dependent occupation and current fractions; in other words, when indeed A occurs, that is, when /It = /i 
and jT = j, then 

^ C+{xt)dt = T J fi{x)C+{x)dx = TTG{n) (30) 
J CQ{xt,dxt) = j dx{G -VU){x) ■ j 6{xt ~ x) o dxt 

= -^[G-jdx (31) 



As a consequence simplifies: the density between the path-space distributions comes out of the path-integral and 
since by construction, for T oo, 1, we get (|27p in the form 

= TFi^l) -Tg{ii) + '^ j{G^ F)-j dx (32) 

That is to be read as the sum of an excess instantaneous traffic given density alone, and an excess work (or, 
equivalently, entropy flux) given stationary current j alone. The excess is being understood in the sense of the above 
modification, with G producing the current j, i.e., j = XP^{G — VJ7) — D^fi, see or 

f3G = {fiD)-\] + \7\og{tie^'') (33) 

Since, by (|23p. excess entropy production equals excess traffic (for our diffusion processes), we can equivalently write 
([5^ as 

= ^^^^^^ + f/(G-^)-,d. (34) 

Upon substituting ([55)) into ([M)) and using p4p . the fluctuation functional can also be written in the following 
explicitly positive form: 



H^^,J) = 1 +2l{^^D) {j-]t,)-jAx 

= \ fij-J^)-(pDr\3-J^)dx 



(35) 



(On the assumption V • = 0; remember that = oo otherwise.) This last formula resembles the Gaussian-like 

expressions for the current distribution, typical for hydrodynamic fluctuations of the diffusion- type. Such expressions 
are omnipresent in the works of e.g. 0, Although the quadratic integrand in (|35p resembles the (generalized 
Onsager-Machlup) Lagrangian for macroscopic fluctuations in the hydrodynamic limit, we have no spatial/temporal 
rescaling here. We have started from a mesoscopic system as described by a diffusion equation and the only large 
parameter is the time span T. What we stress here with respect to other work on dynamical large deviations, 
becomes visible from the formulas (|32p and p4p and has as such nothing to do with hydrodynamic rescaling or with 
macrostatistics. It concerns the thermodynamic interpretation of the fluctuation functional /(/Lt, j) for our mesoscopic 
system: how it is shaped from quantities like traffic, work and entropy production, and providing full account of the 
steady dynamical fluctuations in both the time-symmetric and the time-antisymmetric sectors. 



C. Steady state fluctuation symmetry 

As proof of internal consistency we check the fluctuation symmetry, cf . @, . 
As a corollary of ((55|) . one has 



lifi, -j) - /(m, j) - y J • i^LD)-'j^dx = Sij) (36) 

where 

SU) j F-jdx (37) 
= — y F{xt) o Axt (when jr = j) 

is f3 times the time-averaged power of the nongradient forces under the condition that jx ~ j is the time-averaged 
current. Due to the stationary condition (V • j = 0) it also coincides with the entropy flux. 



D. Result 2: canonical structure of the joint fluctuations 



A useful feature of equilibrium statistical thermodynamics is that the fluctuation functions can in a simple and 
general way be given as differences of thermodynamic potentials that generate the relevant order parameters. As 



suggested already by formula p2p . a similar structure also emerges for the dynamical fluctuations. In order to make 
that manifest we should follow the exact dependence on the driving F, and we therefore indicate that dependence 
here explicitly. In that way, we have an immediate rewriting of (|32p : 

/^^(/i, j) = /o(/^, j) - ^ + rF{^l) - (38) 

where loifJ-ij) is the fluctuation functional for the reference equilibrium process with no driving {F ~ 0), traffic was 
defined in (pij) and the entropy flux follows from (j37p . The canonical structure of that fluctuation functional arises 
because the last difference, the excess traffic 

H{fi,F) = 2[TF{^i)-%{pi)] 
is a potential for the current in the sense that 

Its Legendre transform is 

G(^,j)-sup [SFU)-Hi^l,F)] 

F 

with associated 

^(^,jJ)=/3F(x) (40) 
so that force and current {F, j) make a canonical pair. It can be checked immediately that 

G(M,J) = ^ l{j~f,)-{p^Dr\j~]l)dx (41) 



2 . 

which coincides with 2/o(/i, j) whenever V • j = 0. Using that extended functional to replace Io{iJ,,j) in (|38p. we get 

/f(m,j) = ^[G{fi,j)-SFij)+H{pi,F)] (42) 

for all densities fi and stationary currents j. 
We can still rewrite in the form 

4/j.(m, j) = sup{25j.(j) - (7f(m)} - 2SF{j) + <JF{f^) (43) 

F 



fully in terms of entropic quantities, due to 
A similar canonical structure, cf. ([5^ and (|in|) . has been established already before in the framework of jump processes 
at least on a sufficiently fine-grained scale of description, see [l^. One can therefore conclude that the structure of 
(j42p is also canonical in the sense of being generally valid for a very large class of dynamics. 



E. Small fluctuations 



We look here at the Gaussian approximation in a dynamics far from equilibrium. Later we will also make the 
driving F small, to be close to equilibrium. 

As is clear from ([35]), current and occupations are coupled. It is because of this coupling that contractions of I{p,,j) 
to and to become rather complicated. Even for small fluctuations this coupling remains: take fi ~ p{l + e/ii) 
and j = jp + eji, with e a small parameter. Because j — jf^ is then 0(e), the fluctuation functional is 0{e^): 



(44) 



The last term in this approximation gives the couphng between occupation and current fluctuations. It is proportional 
to the stationary current, which is non-zero away from equilibrium. It is only when we take a dynamics close to 
equilibrium, i.e. F = eFi, that the fluctuations decouple. In this approximation we have that jp = 0(e), and thus, 
near equilibrium, 

^(m, j) = J y dx[ji • {pD)-^j, + Va^i • pDV^Ji^] + 0{e^) (45) 
with, to leading order, a complete decoupling between the time-symmetric and the time-antisymmetric sectors. 



V. CONTRACTIONS 



Now that we have a fluctuation functional for both symmetric and antisymmetric variables, we can compute the 
statistics of empirical averages of arbitrary physical quantities. In particular, we can try to find the fluctuation 
functionals for density /(/x) and and for current separately. The technique to do that is called contraction as we 
go so to speak to a more contracted description. To start we look at the fluctuations of the density (alone). 



A. Result 3: occupation statistics 

Look back at the definition (j25p . As /(/i) = infj I(fi,j), we have to compute the minimizing current j for any given 
density fi. Since the minimization is constrained via the stationary condition V • j = 0, we get the equation 

J = X/i [J^ - V • ([/ + i')] - DV^l (46) 

where tp is a, Lagrange multiplier (function of x). Not surprisingly, we see that the minimizer is the stationary current 
for a modified dynamics that makes fi stationary. This modified dynamics is achieved here by replacing the imposed 
potential U with a modified one, called V = U + tjj (for some ip). We therefore call the minimizing current in (j46p 
j^, and the fluctuation functional becomes: 

m 01 - J,) ■ il^DyHj^ - J,) dx (47) 

For some explicit examples of solutions to see furtherdown, in equations 

The fluctuation functional obtains other equivalent forms by following a road backward from the one that led 
us before to (^5]) : 

/(^) = T(M) - Mm) = ^^^^^^ (48) 

where the second equality follows again from (|23p . In this way we have recognized the excess traffic (or here also: 
the excess entropy production) as governing the large time statistics of the occupation times. That constitutes our 
third main result: the fluctuation functional I{fJ.) exactly equals one quarter of a difference in entropy production 
rates when having density /i, these rates being computed respectively for the original dynamics and for a modified 
dynamics that makes n stationary. 



In formula; ([47]) and ([48|) the potential V has to be determined from fj, by solving the inverse stationary problem (|46p . 
We now give two classes of examples where this V and hence can be made explicit. 

Diffusion on the circle. For the one-dimensional example of Section [II Al the inverse stationary problem (j46p allows 
for an explicit solution. The current jj^ is immediately read off the formula ([7]) , 



3l - .1 , . , / FAX (49) 







and the potential Y obtains the form 



1 r 7^ 



Vix)^--lo,,ix) + J^ {F^^)Ay (50) 



which is a nonlocal functional of the given density /i. The fluctuation functional is explicitly given as 



AI{^,)=ai^,)^ (51) 
Jo {f^D)-^dx 

for n^O. 

Observe that if /i = on some open set A then the rate function equals /(/i) = (T(/i)/4. (That follows also from the 
equilibrium form (j54[) below as the circle gets effectively cut and the dynamics mimics a detailed balance one.) The 
infimum of I{p) over all densities fi that vanish on A then gives the escape rate from the complement A"^ — [0, 1) \ A. 
As a simple example, assume that U ~ and let F and D be some constants. In this case the entropy production 
(fU]) reads 

CT(/i) = P^DF^ + D r ^dx (52) 
Jo M 

To compute the escape rate from A'^ (or, entrance rate to A) we must take the infimum of ([5^ over all /i that vanish on 

A. Setting A = (0, S) for some < i5 < 1, that infimum is reached for the density fi*{x) = sin^(^^^3^), x G [S, 1], 
and the escape rate is 

tt'^D B'^DF'^ 

inf /(m) - I{n - j--^ + '—— (53) 
mU=o (1 - oy 4 

Even in equilibrium {F = 0) the result is meaningful as it relates the diffusion constant to an escape rate. In the 
context of dynamical systems, the analysis of the escape rates and of their link to linear transport coefficients was 
initiated by Dorfman and Gaspard, sec [26. . 27\ and references therein. 

Close-to- equilibrium dynamics. Let us start from what can be said in general for equilibrium diffusions. If _F = 
then equation (|46p has the solution V = —13~^ log ^i, and the corresponding current and the entropy production 
(Jv{l^) are both zero. As a result, 

= ^ (54) 

This exact relation between the equilibrium dynamical fluctuations and the entropy production is solely true for 
diffusion processes. In contrast, for jump processes gives only the leading term in an expansion of around 
the equilibrium density p oc e~^^ , and the relation (|54p obtains corrections when beyond small fluctuations; see 
for details. 

When close to equilibrium, the computation of /(/i) by contraction is easy: we see from that the second term 
on its right-hand side is just /(/i). 

In the same approximation of small fluctuations and close-to-equilibrium, the entropy production becomes: 

ai^l) = J dx[jp ■ {pD)-\jp + e^VMi • pDV^il] + 0(6^) (55) 

and thus we get 

^^^^^ + 0(e3) (56) 

This reveals to be a special case of a general result, according to which the entropy production governs the 
occupational statistics in the linear irreversible regime. It provides a fluctuation-based explanation for the minimum 
entropy production principle introduced by Prigogine to characterize stationarity via an (approximate) variational 
principle, (28j : the stationary state has minimal (not necessarily zero) entropy production. 



B. Current statistics 



The contraction to the current j is also possible. However, up to special examples, there is no explicit solution to 
the associated variational problem and for general models one has to resort to a perturbative or numerical analysis. 
In fact, often the calculation starting from the generating function of the current appears more practical than to do 



the contraction starting from see (ol. [Tl|. 

We restrict us here to giving the resuh for a constant drift on the circle and to smaU fluctuations around equihbrium. 

Constantly driven diffusion on the circle. Again we take U = and F, x constants. In this case, from (|35|1 the 
joined fluctuation functional reads: 

^(m, j) = ^ / ^ (j - PDF^i -D^i'f dx (57) 
and for all j, the inflmum over n is reached at the uniform distribution, so that 

,0) . iq^ ,58, 

and hence we see that here the current fluctuations are Gaussian. 

Close-to-equilibrium. We have the analogue of the minimum entropy production principle. The starting point is 
again from which we extract the current fluctuations: 



lU) = 7/0 - Jp) ■ {PD)-\J - Jp) Ax + 0(63) 

Y (59) 
^-[&{jp) + &U)-2S{j)\+0{e^) 



with 2>{j) ~ J j ■ {pD)^^j da; sometimes called the Onsager dissipation function, and S{j) is the entropy flux defined 
in In particular, this leads to a variational characterization of the steady current jp which can be written as the 
following maximum entropy production principle: the jp maximizes the entropy flux S{j) under the two stationary 
constraints 

(l)V-j = 0, (2)^(j)=5(j) (60) 

The second condition is indeed satisfied at j = jp (note also that p can with no harm in this order be replaced by 
the equilibrium density pf=o = e~^^ /Z.) Such a variational principle, known as a maximum entropy production 
principle, is often used in applications and apparently even beyond the linear irreversible regime. As is however 
clear from (j59p from our dynamical fluctuation theory, the validity of the maximum entropy principle is restricted to 
close-to-equilibrium. Beyond that regime, we must refer to contractions from (j32p . (|42p or even from (j44p for generally 
valid expressions with a general thermodynamic meaning. 
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APPENDIX A: STOCHASTIC INTEGRALS 

We collect here some necessary albeit formal manipulations from stochastic calculus. We refer to [1^ and for 
further systematics. 

Physical quantities, such as work and heat dissipation, depend on the speciflc trajectory (or path) that the system 
covers during its evolution. But for diffusions these paths are not differentiable. Therefore integrals like Jg f{xt)Axt 
cannot be deflned in the usual way. It turns out that these integrals, called stochastic integrals, can be interpreted in 
different ways. Most common are the Ito and Stratonovich interpretations, see e.g. [29l |. 

ltd integral. For the Ito interpretation, the integral domain [0, T] is split up in a set of discrete points = io < 
ti < . . . < tn = T, with Atj = tj — tj-i, such that At = max-,- Atj for n — > 00. It is important to note that for 
diffusions ([5]) we have that for Atj 0, 



(Al) 



where a and /3 denote the components of the vectors and the diffusion matrix. The stochastic integral is then computed 
as 

/ f{xt)dxt^ hm V/(a;t^._J-(xt^-Xt^._J (A2) 

JO 

Note that the function / can be either scalar or vector, and it is evaluated at the left endpoint of the intervals [tj-i, tj\. 

For the Ito integral one cannot use the normal rules of integration. Instead it is easily checked from (|Aip that for 
any function g, 

[ Vg{xt) ■ dxt = g{xT) - g{xo) [ (DV ■ Vg){xt) dt (A3) 
Jo Jo 

(The symbol ■ stands for the scalar product.) 

Stratonovich integral. The Stratonovich interpretation differs from the Ito interpretation only in the points of 
evaluation of the function /. In this case / is evaluated in the midpoints of the time intervals: 

n X -\r X 

^ (A4) 

where the symbol o is commonly added as a notation to distinguish between Ito and Stratonovich interpretations. 
For our analysis it is important to observe that the Stratonovich integral is time-antisymmetric, and also that 

T 

Vgixt) o dxt = g{xT) - gixo) (A5) 



Relation between ltd and Stratonovich. For a scalar /, it is easily found from (jAl|) that 



f{xt)odxt= f{xt)dxt+ {DVf){xt)dt (A6) 
Jo Jo 

(and analogously for vectors). Using that xt solves the Ito-stochastic equation (jS]), that Stratonovich integral further 
explicitly equals to 

/ f{xt)odxt^ j [fx{F-VU) + \/ ■{fD)]{xt)dt+ [ f{xt)^2D{xt)dBt (A7) 
Jo Jo Jo 

Observe that the last Ito-integral has mean zero since the integrand evaluated at each mesh point Xt _-^ and the 
increment Bt- — Bt-_^ of the Brownian motion are mutually independent, and the latter has zero mean. Hence, the 
mean value of the Stratonovich integral (|A6|) is 

( r f{xt) o dxt) = { CuxiF - VC/) + V • {fD)][xt) dt) 

^Jo ' !-'■« ^JO ' f^o 

^ dt J ^it [fxiF -VU) + \7. ifD)] dx (A8) 
= At J f [x{F - yU)^lt - DVfit] dx 



The first equality is obtained by taking the mean of (IA7I1 . The second equality uses that /i^, defined by the Fokker- 
Planck equation is the evolved density at time t when starting from fiQ at time zero. Therefore Xt is there 
distributed according to fif The final equality is from a partial integration. That proves the equality (jlip . 



APPENDIX B: PATH-DISTRIBUTION OF DIFFUSION PROCESS 



The path of a particle subject to the stochastic equation ([S]) is completely determined by the initial condition and 
the realization of the Brownian motion Bf. Since the latter is a standard Gaussian process, we can integrate out the 
noise to obtain a path-integral representation. We restrict here to a simple heuristic argument which goes as follows. 

To find the probability for the particle to pass through (the neighborhood of) a discrete set of points lu — (xq, to = 
0; xi, ti; . . . ; a;„, i„ — T), and having already in mind the limit n — > oo, we use that each increment Axj = Xj — Xj^i 
is by ([H]) in a one-to-one correspondence with an increment of the Brownian motion, which is {2D{xj-i))~^^'^[Axj — 
[xF — xVC/ -|- V • D) Atj]. Since the latter increments are independent standard Gaussian variables, we immediately 
obtain the discrete path-distribution (or its density with respect to the [n + 1)— product of flat distributions): 

n 

dP^„(w) - n{xo)5xa JJ[47rD(xj_i) At^]-^ 

X exp{-^[Ax, - ixF - xVC/ + V ■ D){xj-,) At^] 

■ D-\xj^i)[Axj - (xF -yu + V ■ D){x^^i)At^] | 5xj (Bl) 

where the initial point xq is sampled from an initial distribution fiQ. A useful representation employs the reference 
process 

dxt = V ■ D{xt)dt + yj2D{xt)dBt (B2) 

with the initial point sampled from the flat distribution; note the latter (unnormalizable in general) distribution is 
stationary. Its path-distribution P° is simply obtained from (|B1[) by putting F = VU = and /i = 1. Comparing 
both path-distributions and passing to the limit n —>■ oo, At — > 0, we obtain a Girsanov-formula for diffusions: 



r 3 f 

^(u;)=/x(xo)exp|-^ [{F - VU) ■ xiF ~ yU)]{xt) dt 



dp' 



where an Ito integral comes out by construction. Using the relation (jA6|) to replace that integral with a Stratonovich 
one, we finally arrive at formulas (jl5p - (|17p . A rigorous derivation can e.g. be found in [30| . 
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